EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-PH-EP/2005-00 1 
18 January 2005 
revised^ 22 March 2005 



Measurement of the Spin Structure of the Deuteron in the 

DIS Region 

The COMPASS Collaboration 



Abstract 



We present a new measurement of the longitudinal spin asymmetry Af and the spin- 
dependent structure function gf of the deuteron in the range 1 GeV 2 < Q 2 < 100 GeV 2 
and 0.004 < x < 0.7. The data were obtained by the COMPASS experiment at CERN 
using a 160 GeV polarised muon beam and a large polarised 6 LiD target. The results are in 
agreement with those from previous experiments and improve considerably the statistical 
accuracy in the region 0.004 < x < 0.03. 
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Since the surprising result obtained for the spin structure function of the proton by the EMC Q, 
the determination of the longitudinal spin structure of the proton and the neutron has remained 
one of the important issues in particle physics [2|. The spin structure functions are used to test 
the Bjorken sum rule and to determine quark and gluon polarisations from the QCD evolution 
equations [3]. They are also used as constraints in the derivation of the polarisation of quarks of 
different flavour from semi-inclusive asymmetries 013. 

Here we report on the first results from the COMPASS experiment at CERN on the 
deuteron spin asymmetry Af and the spin-dependent structure function gf in the deep inelastic 
scattering (DIS) region, covering the range 1 GeV 2 to 100 GeV 2 in the photon virtuality Q 2 and 
0.004 to 0.7 in the Bjorken scaling variable x. 

The COMPASS spectrometer is located in the same muon beam line as the former SMC 
experiment and covers a similar kinematic region for inclusive reactions. However, it uses a 
higher intensity muon beam of 160 GeV, a longitudinally or transversely polarised target made 
of 6 LiD, and a new two-stage spectrometer. A general description of the experiment has been 
presented in Ref. [6| and only the most relevant elements for the present analysis will be men- 
tioned below. The data in the longitudinal configuration taken in 2002 and 2003 correspond to 
luminosities of about 600 pb _1 and 900 pb _1 , respectively. 

The experiment was performed at the M2 muon beam line of the CERN SPS. The muons 
originate from the decay of n and K mesons produced by the 400 GeV proton beam on a 
primary beryllium target. The fi + intensity is 2 • 10 s per spill of 4.8 s with a cycle time of 
16.8 s. The beam profile presents a Gaussian core and a large non-Gaussian tail due to halo 
muons. The beam has a nominal energy of 160 GeV and is focused at the target centre, with 
a spread of 7 mm (r.m.s.) and a momentum spread of a p /p = 0.05 for the Gaussian core. The 
momentum of each muon is measured upstream of the experimental area in a beam momentum 
station consisting of five (four in the year 2002) planes of scintillator strips with a dipole magnet 
in between. The precision of the momentum determination is typically Ap/p = 0.003. The 
incoming muon direction and position are measured by small scintillating fibre hodoscopes 
and silicon microstrip detectors J7J0D. The space resolution is about 0.12 mm for the fibres 
and 0.015 mm for the microstrips, and the direction of the incoming muon is measured with a 
precision of 30 /xrad. 

The polarisation Pb of the beam muons was determined by a Monte Carlo program mod- 
elling in detail the phase space of the parent hadrons and decay muons, as well as their prop- 
agation through the beam transport system 0. Within a precision of about 0.04 the calculated 
values are consistent with the polarisation measurements performed by the SMC at 100 and 
190 GeV [ 10 1 . For the present experiment the model gives a polarisation of the muon varying 
with its energy from —0.57 at 140 GeV to —0.86 at 180 GeV with a mean value of —0.76. 

The target is located inside the solenoid magnet previously used by the SMC experiment 
ifTTl . which provides a field of 2.5 T along the beam direction. The magnet aperture seen from 
the upstream end of the target is ±70 mrad. The target consists of two cells, each 60 cm long and 
3 cm in diameter, separated by 10 cm. They are filled with 6 LiD which is used as deuteron target 
material and longitudinally polarised with dynamic nuclear polarisation (DNP) lfT2l . The two 
cells are polarised in opposite directions so that data from both spin directions are recorded at 
the same time. The polarisation is measured by NMR coils with a relative precision of about 5% 
|[T3l . The typical polarisation values obtained after a build-up time of about 5 days are +0.53 
and —0.50. The spin directions in the two target cells are reversed every 8 hours by rotating the 
magnetic field direction. In this way, fluxes and acceptances cancel out in the calculation of spin 
asymmetries, provided that the ratio of acceptances remains unchanged after spin reversal. In 
order to minimise possible acceptance effects related to the orientation of the solenoid field, the 
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Figure 1: Layout of the COMPASS spectrometer used in 2003. The configuration was identical 
in 2002 except for Electromagnetic Calorimeter 2 which was not included in the read-out. The 
thin vertical lines represent the tracking detectors. 



sign of the polarisation in each target cell is also reversed several times per year by changing 
the DNP microwave frequencies. 

The COMPASS spectrometer (Fig.[U) is designed to reconstruct the scattered muons and 
the produced hadrons in wide momentum and angular ranges. It is divided in two stages asso- 
ciated with two dipole magnets, SMI and SM2. The first one is a large- aperture magnet, with a 
field integral of 1 Tm along the beam line, which accepts charged particles of momenta larger 
than 0.4 GeV. The second magnet, SM2, has a field integral of 4.4 Tm and accepts particles 
of momenta larger than 4 GeV. Different types of tracking detectors are used to cope with the 
rapid increase of the particle rate from the outside to the central beam region. The beam region 
downstream of the target is covered by scintillating fibre detectors 0, the region near to the 
beam by micromesh gaseous chambers iTPfl and gas electron multiplier chambers [ 15 1. The in- 
termediate region, further away from the beam line, is covered by drift chambers and multiwire 
proportional chambers. Large-angle tracking is mainly provided by straw detectors |[T6ll and by 
large drift chambers. The identification of muons is based on the fact that they are observed 
behind hadron absorbers. Two 'muon wall' detectors are used: the first one, located in front 
of SM2, consists of two stations of Iarocci-type chambers with an iron layer in between and 
detects muons outside the aperture of SM2; the second one, installed at the end of the spectrom- 
eter, is composed of drift tubes and detects the muons which passed through SM2. Hadrons are 
detected by two large iron-scintillator sampling calorimeters, installed in front of the absorbers 
and shielded to avoid electromagnetic contamination. 

The data recording system is activated by a combination of signals indicating the presence 
of a scattered muon at a given angle or in a given energy range. In most DIS events (Q 2 > 
1 GeV 2 ), the scattered muon is identified by coincidence signals in the trigger hodoscopes, 
that define its direction behind SM2. Several veto counters installed upstream of the target are 
used to avoid triggers due to halo muons. In addition to this inclusive trigger mode, which 
was commonly used in the previous CERN muon experiments, several semi-inclusive triggers 
select events fulfilling requirements based on the muon energy loss and on the presence of a 
hadron signal in the calorimeters ifTVl . Calorimeter signals due to halo muons are rejected by 
requiring the presence of at least one cluster with an energy deposit exceeding three times the 
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Figure 2: Fraction of inclusive, semi-inclusive, and calorimetric triggers in the final data sample 
(2002 and 2003) as a function of x (left) and Q 2 (right). Events are counted with the weight 
they carry in the asymmetry calculation [Eq. (7)]. 



average value expected for a muon. This condition provides a trigger efficiency of more than 
80% for events with total hadronic energy E had > 30 GeV. In a part of the 2003 data taking, 
the acceptance was further extended towards high Q 2 values by the addition of a standalone 
calorimetric trigger in which no condition is set for the scattered muon but an energy deposit 
in the hadron calorimeter exceeding 9 times the typical muon response is required. The semi- 
inclusive and calorimetric triggers thus select a sample of hadronic events which are analysed 
in parallel with the inclusive sample. The relative contributions of the different trigger types 
are shown as a function of x and Q 2 in Fig. |2j The fraction of inclusive triggers, where the 
selection criteria refer only to the scattered muon, varies from 60% to 75% over the range of x 
(events satisfying simultaneously inclusive and non-inclusive trigger conditions are counted as 
inclusive). The semi-inclusive triggers account for about 40% of the data at low x and decrease 
steadily for x > 0.02, while the contribution of the standalone calorimetric trigger starts around 
x = 0.02 and reaches 30% in the highest x bin. 

Larger variations of the different contributions are observed as a function of Q 2 : the in- 
clusive triggers account for 80% of the events at medium Q 2 (3-15 GeV 2 ), while the standalone 
calorimetric trigger becomes dominant for Q 2 > 30 GeV 2 . 

In order to eliminate spurious triggers as well as badly or partially reconstructed events, 
a reconstructed interaction point connected to a beam muon and to a scattered muon is required 
for all events. In addition, the presence of a hadron track at the interaction point is required for 
the semi-inclusive and standalone calorimetric triggers. The track reconstruction efficiency was 
found to be about 95% for scattered muons and for high-energy hadrons (E > 30 GeV) that 
were generated in a Monte Carlo simulation, tracked through the spectrometer, and analysed 
in the same way as the data. The direction of tracks reconstructed at the interaction point is 
determined with a precision better than 0.2 mrad and the momentum resolution for scattered 
muons is about 0.5%. 

As the COMPASS trigger setup is predominantly intended for the study of quasi-real 
photon interactions, DIS events represent only a small fraction of the data sample. The combi- 
nation of cuts on the photon virtuality (Q 2 > 1 GeV 2 ), the fraction of energy carried away by 
the virtual photon (0.1 < y < 0.9), and the requirement that the interaction take place within 
one of the target cells results in a reduction factor of about 20. In addition, the incoming muon 
momentum is required to be in the interval 140 GeV < < 180 GeV and, in order to equalise 
fluxes seen by the two target cells, its trajectory is required to cross entirely both target cells. For 
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Figure 3: The dilution factor / of the LiD polarised target as a function of x for inclusive and 
hadronic events. The dilution due to radiative effects on the deuteron is included. The values 
quoted for each x bin are averaged over the kinematic range of the corresponding triggers. 



consistency, in events triggered by hodoscope signals, it is also verified that the reconstructed 
scattered muon hits the hodoscopes that have generated the event trigger. The resulting sample 
amounts to about 34 • 10 6 events with a fraction of 71% of the data collected in 2003. 

The cross-section asymmetry A d = (<t^ — a^)/(a^ + cr^), for antiparallel (f|) and 
parallel (f f ) spins of the incoming muon and the target deuteron, is related to the virtual-photon 
deuteron asymmetries Af and A d by 

A d = D(Ai + v A d 2 ), (1) 

where the factors 77 and D depend on the event kinematics. The virtual-photon depolarisation 
factor 

D y{2 ~ y) (2) 

depends in addition on the unpolarised structure function R = o~l/(j t . The longitudinal virtual- 
photon deuteron asymmetry is defined as 

At = {al - a T 2 )/{2a T ) , (3) 

where crj is the virtual-photon-deuteron absorption cross-section for total spin projection J in 
the photon direction, and a T = (1/3) (<Jq + + a%) is the total transverse photo-absorption 
cross-section. The transverse asymmetry A d has been accurately measured ifTHI and was found 
to be small. Since the kinematic factor 77 = ^-y) VQ^/E^ is also small in the COMPASS 
kinematic range, the second term in Eq. (1) can be neglected, so that 

Af ~ A d /D , (4) 

and the longitudinal spin structure function is given by 

pd 

91 = 2x(l + R) Adl ' (5) 

where F d is the deuteron spin-independent structure function. The number of events Ni col- 
lected from a given target cell in a given time interval is related to the spin-independent cross- 
section a and to the asymmetry Af by 

Ni = airmail + P B P T fDAf) , (6) 
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Figure 4: The asymmetry Af(x) as measured in COMPASS and previous results from SMC 
ESI . HERMES 0, SLAC E143 $M and E155 \21\ at Q 2 > 1 GeV 2 . The SLAC values of 
gi/Fi have been converted to A\ and the El 55 data corresponding to the same x have been 
averaged over Q 2 . Only statistical errors are shown with the data points. The shaded areas show 
the size of the COMPASS systematic errors. 



where Pb and Pt are the beam and target polarisations, fa the incoming muon flux, aj the 
acceptance for the target cell, rii the corresponding number of target nucleons, and / the target 
dilution factor. For a 6 LiD target the dilution is naively expected to be of the order of 50% 
because 6 Li can be described as an 4 He core and a deuteron [ 19 1. The dilution factor / is given 
by the ratio of the absorption cross-sections on the deuteron to that of all nuclei entering the 
target cells. It includes a correction for the relative polarisation of deuterons bound in 6 Li with 
respect to free deuterons. It also includes the dilution due to radiative events on the deuteron, 
which is taken into account by the ratio of the one-photon exchange cross-section to the total 
cross-section p = <J d l<j\ ot lEUll . The values of / are shown in Fig. as a function of x for 
inclusive and hadronic events. The large difference observed at low x results from the factor p 
which is much smaller in the inclusive case because radiative effects in elastic scattering largely 
contribute in the denominator. The dilution factors also differ slightly at high x because the 
inclusive and standalone calorimetric triggers cover different ranges of Q 2 as shown in Fig. El 

The asymmetry is extracted from data sets taken before and after a reversal of the tar- 
get spin directions. The four relations of Eq. (6), corresponding to the two cells (u and d) 
and the two spin orientations (1 and 2), lead to a second-order equation in A\ for the ratio 
(N Ut iN dj2 ) I (N dj iN u>2 ). Fluxes and acceptances cancel out in this equation if the ratio of accep- 
tances for the two cells is the same before and after the reversal iPTl. In order to minimise the 
statistical error, all quantities used in the asymmetry calculation are evaluated event by event 
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Table 1: Values of Af and gf with their statistical and systematical errors as a function of x 
with the corresponding average values of Q 2 and y. 
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with the weight factor 

w = P B fD. (7) 

The polarisation of the beam muon, P B , is obtained from a parametrisation as a function of the 
beam momentum. The factors / and D are calculated from the kinematic variables with the 
value of R taken from the NMC |22| or the SLAC parametrisation [23| for x below or above 
0.12, respectively. The target polarisation is not included in the event weight [Eq. (7)] because 
it may vary in time and generate false asymmetries. An average Pt (~ 0.5) is used for each 
target cell and each spin orientation. 

Inclusive and hadronic events are analysed separately with the corresponding value of the 
dilution factor. The additive radiative correction to the asymmetry [21 ] has also been calculated 
separately ll24ll using an input parametrisation of Af fitted to the present data. The values ob- 
tained for inclusive and hadronic events differ by 0.0003 in the lowest x-intervals and become 
nearly equal at higher x. These additive corrections are negligible at low x and reach a maxi- 
mum value of 0.008 at high x. 

The asymmetries obtained for hadronic events are statistically compatible with the inclusive 
ones and their differences do not show any hint of a systematic dependence on x. This obser- 
vation agrees with the Monte Carlo study of Ref. ll25l which also shows that the selection of 
hadronic events has no sizeable effect on the evaluation of Ai for interactions on a deuteron 
target within the kinematic range and the hadron acceptance of the present experiment. 

The final values of Af are obtained by merging the inclusive and hadronic sets weighted 
according to their statistical errors. They are listed in Table 1 with the corresponding statistical 
and systematical errors and shown in Fig. |4|in comparison with those obtained by the SMC 
ESI , by E143 m and E155 \27\ at SLAC, and by HERMES [5|. Good agreement is observed 
over the full range of x. For the four points with x < 0.03, our results reduce the statistical 
errors of previous measurements by a factor of about 2.5. 

Figure |5] shows the values of Af as a function of Q 2 for each interval of x. The results 
of fits to a constant in each interval of x are shown by the solid lines. They yield an average 
X 2 -probability of about 0.5 and do not indicate any Q 2 dependence. Some dependence of Af on 
Q 2 is expected from perturbative QCD because the Q 2 evolutions of spin dependent and spin 
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Figure 5: Values of Af as a function of Q 2 in intervals of x. The solid lines are the results of fits 
to a constant; the dashed lines show the Q 2 dependence predicted by perturbative QCD. 



independent structure functions are different. However previous experiments ll25ll have shown 
that the two Q 2 evolutions largely cancel out so that the values of Af at fixed x become nearly 
independent of Q 2 . The Q 2 dependence predicted by the SMC fit of Ref. [28| is shown by the 
dashed lines in Fig.|5]and describes the data equally well. 

The systematic error on Af contains an overall scale uncertainty of 6.5% due to the un- 
certainties on Pb and Pt. The error on the dilution factor /, which takes into account the uncer- 
tainty on the target composition and the uncertainty on the corresponding cross-section ratios, is 
of the order of 6% over the full range of x. The uncertainty on the parametrisation of R affects 
the depolarisation factor D [Eq. (2)] by 4-5%. The neglect of the A 2 term mainly affects the 
highest x interval where its contribution is estimated to be < 0.005. The error on the radiative 
corrections to the asymmetry is estimated by varying the input parametrisation of Af(x) within 
the statistical error of the present data. The effect of event migration to neighbouring x bins, 
resulting from the smearing of kinematic variables due to the finite resolution of the spectrom- 
eter and to the radiative effects, was evaluated by a Monte Carlo simulation and found to be 
negligible. Potential false experimental asymmetries were searched for by modifying the selec- 
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Figure 6: Values of x gf(x) vs. x. The COMPASS points are given at the (Q 2 ) of each interval 
of x. The SMC points ll25ll were evolved to the Q 2 of the corresponding COMPASS point and 
are slightly shifted to larger x for clarity. Only statistical errors are shown with the data points. 
The upper and lower bands show the COMPASS and SMC systematic errors, respectively. 

tion of data sets used for the asymmetry calculation. The grouping of data into configurations 
with opposite target-polarisation was varied from large samples, covering at most two weeks of 
data taking, into about 100 small samples, taken in time intervals of the order of 16 hours. A 
statistical test was performed on the distributions of the asymmetries extracted from these small 
samples. In every interval of x they were found to be normally distributed, with a standard de- 
viation a compatible with the one derived from the statistical errors (a sta t)- Time-dependent 
effects which would lead to a broadening of these distributions were thus not observed. Since 
the spread of the observed cr's is about 0.05, we take 1.1 a sta t as upper limit for a and obtain 
for each x bin a conservative upper bound of the systematic error arising from time-dependent 
effects 



Asymmetries for configurations where spin effects cancel out were calculated to check the can- 
cellation of fluxes and acceptances. They were found compatible with zero within their statis- 
tical errors. The comparison of asymmetries obtained from different parts of the spectrometer 
did not show any systematic effect. Asymmetries obtained with different settings of the DNP 
microwave frequency were compared in order to test possible effects related to the orientation 
of the target field. No sizeable effect was observed. 

The values of gf(x, Q 2 ) quoted in the last column of Table 1 were obtained from Eq. (5), 
with the F$ parametrisation of Ref . [ 25 1 and the parametrisation of R already used in the calcu- 
lation of the depolarisation factor. The systematic errors on gf contain an additional contribution 
due to the uncertainty on the parametrisation of F£. The error due to the uncertainty on R is re- 
duced by a partial cancellation between the R dependence of the depolarisation factor [Eq. (2)] 
and the factor (1 + R) in Eq. (5). Our values of gf are shown in Fig. |5]in comparison with 



Csyst < 0.5 (Jstaf 



(8) 



8 



the SMC results [25| which cover the same Q 2 range and were evolved to the same Q 2 values. 
Their improved accuracy provides a better evaluation of gi at low x: integrating the values of 
gf(x) shown in Fig.|6]over the range 0.004 < x < 0.03, we obtain (—0.3 ± 1.0) ■ 10 -3 and 
(-5.3 ± 2.3) ■ lO" 3 for COMPASS and SMC data, respectively. For x < 0.03 the COMPASS 
results are consistent with zero and do not show the tendency of the SMC data of negative gf 
values. 

In combination with the accurate SLAC and HERMES data at larger x, our new results 
will improve the extrapolation of gf towards x = 0. However, taken alone, they do not provide 
a more accurate evaluation of the first moment Tf because of the relatively large errors at high 
x resulting from the late implementation of the calorimetric trigger in the present data. These 
errors will be reduced for the 2004 data where the calorimetric trigger was used during the full 
data-taking period. 

In conclusion, a new evaluation of the longitudinal spin asymmetry and the spin structure 
function of the deuteron in the DIS region (Q 2 > 1 GeV 2 ) was performed by the COMPASS 
experiment at CERN. The data cover nearly the same range of x as the former SMC experiment, 
0.004 < x < 0.7. The results are in agreement with previous experiments over the full range of 
x and significantly improve the statistical accuracy in the region x < 0.03. 
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